Effective size of a trapped atomic Bose gas 
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We investigate the temperature-dependent effective size of a trapped interacting atomic Bose gas 
within a mean field theory approximation. The sudden shrinking of the average length, as observed 
in an earlier experiment by Wang et al. [Chin. Phys. Lett. 20, 799 (2003)], is shown to be a good 
indication for Bose- Einstein condensation (BEC). Our study also supports the use of the average 
width of a trapped Bose gas for a nondestructive calibration of its temperature. 
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I. INTRODUCTION 

It has been a decade since the first realization of 
Bose-Einstein condensation (BEC) in alkali atomic gases 
[E IE 0] • Many properties of atomic quantum gases have 
been extensively studied and understood [E IE 13 • A 
consistent theme has been the general agreement be- 
tween theoretical predictions and experimental obser- 
vations. For example, at low temperatures when the 
non-condensed thermal component is negligible, Gross- 
Pitaevskii equation has been successfully used to model 
the condensate mean field, which provides a firm theoret- 
ical framework for the description of condensate den sity 
distribution and collective excitations [E |E 0, 0, 0, 
HE HE HE HE HE Hll At nonzero temperatures, several 
mean field theory extensions such as the Hartree-Fock- 
Bogoliubov or the Hartree-Fock-Popov approaches gen- 
erally become sufficient, as evidenced by the calculated 
temperature-dependent condensate fraction, ex pan sion 
ene rgy , and number fluctuation, etc. Il9l I20L l2lll22l I23I 

EE lH EE IH El Hi HE EE Hi EE illsElir 

In this paper we report our investigations of the 
temperature-dependent average width of a trapped in- 
teracting Bose gas within the mean field Hartree-Fock 
approximation. We undertake this study in an effort to 
provide some theoretical support for the recent experi- 
ment by Wang et al., which produced the first atomic 
condensate inside mainland China 0. Few earlier the- 
oretical studies have addressed the issue of an averaged 
length/ width of a trapped interacting Bose gas |2^], de- 
spite it being one of the topics studied in the first gener- 
ation of BEC experiments • As expected, our study 
points to an universal feature for all trapped interacting 
Bose gases. We find that the effective, or the average 
size of a repulsively interacting trapped Bose gas drops 
suddenly when condensation occurs, i.e. after the gas 
temperature is lower than the transition temperature Tq- 
As a consequence we support the use of this reduction in 
effective sample size as a good evidence for the onset of 
BEC, distinctively different from the more conventional 
bimodal density profile based on a measurement of the 
expanded condensate density distribution after the con- 
fining trap is switched off. 



This paper is organized as follows. The present sec- 
tion provides a brief introduction. Sections [H] and IIIII 
describe the calculations for the average size of an ideal 
or interacting trapped Bose gas respectively. We analyze 
the energy partitions of a Bose gas in Sec. IIVI Section IVI 
presents a short conclusion which summarizes our results. 

II. SIZE OF AN IDEAL BOSE GAS 

For an ideal Bose gas, the average number of particles 
in a single particle state z) with energy e,; is given by the 
familiar Bose-Einstein distribution 

where f3 = l/(fcsT) with ks denoting the Boltzmann 
constant, g± is the degree of degeneracy for state \i), z = 
exp(/3/x) is the fugacity, and fi is the chemical potential 
determined by the conservation of total number (N) of 
particles Y^tLo Ui ~ ^- 

The statistical properties of a Bose gas, such as its 
specific heat, the condensate fraction, etc, are completely 
determined once the chemical potential is found |19t l2ll 
HE l34l l36j | . For a spherically symmetric harmonic trap 
V cxt (r) — MuV/2, the width for state \i) is easily found 
to be 

where a r = \J 'Ti / \M ~lo) is a characteristic length for the 
harmonic trap. The width of a Bose gas is then obtained 
as 




For large N, we take the usual approximation of changing 
the summation into an integral weighted by the density 
of states and separate out the atoms in the ground state. 
We find that 

<r 2 > ~ Qa 2 )7V + 3a 2 g 4 (z), (4) 



2 



1.5r 



o 1 



0.5 



o A o, + o + <i + * 



0.5 



T/T„ 



1.5 



FIG. 1: (Color online) The temperature-dependent width of a 
Bose gas in a spherically symmetric harmonic trap. The solid, 
dashed, dotted lines denote respectively the case of an ideal 
gas with N = 5 x 10 5 , 1 x 10 5 , and 1 x 10 4 . The correspond- 
ing results for an interacting gas are given by empty circles, 
asterisks, and plus signs. The trap parameters are given in 
Table |3 Note that T c and r c both depend on N. 
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FIG. 2: (Color online) The same as in Fig. 0but for a cigar- 
shaped harmonic trap (the left column) or a pancake-shaped 
harmonic trap (the right column). The total number of atoms 
is assumed iV = 5 x 10 with trap parameters as given in Table 
HI The solid lines are for an ideal gas while the empty circles 
refer to an interacting 87 Rb gas. 



if the temperature satisfies fc^T 3> Jko, where 



N = 



0. 

N - 



T>T C , 
^f)% 3 (z), T<T C , 



with g v (z) = ZZi* k /k u , CM = 9,(1), and T c = 
(hLo/k B )(N/C(3)) 1 / 3 . This expression takes a famil- 
iar form with the first term denoting the squared width 
for the ground state (condensate), while the second term 
for the excited states (thermal component). 

We note there exist two apparent temperature scales 
in the system, Tiu) and k B Tc- When k B Tc S> tuv, as 
normally is the case, our approximation is expected to 
work well. For high temperatures k B T 3> fito, we ignore 
the contribution from the condensate and obtain 



3a„ 



k B T 



T 



(5) 



for T < Tc- On the other hand, we find 



3aiN 



k B T\ g 4 (z) 



huj ) 53 (z) 



ar 2 c 



T 



*r > (6) 



for T > T c . r 2 c = 3a^(4) [7V/C(3)] 4 / 3 denotes the 
squared width at T c and a = 3 4 (z)^(3)/[3 3 (z)^(4)] ~ 1 
has a very weak dependence on T. We see that in a spher- 
ically symmetric harmonic trap, the squared width of an 
ideal Bose gas is proportional to T 4 for T < Tc and to T 
for T > Tc, as clearly shown in Fig. ^ consistent with 
earlier experimental reports that the area of absorption 
image of a Bose gas is proportional to its temperature in 
the absence of a condensate Q, 1^3 . 

For a cylindrically symmetric trap with oj x = u> y = 
u>j_, the temperature dependence of the three widths 
squared is the same as in a spherically symmetric 



trap discussed above. The prefactors become respec- 



tively 



-G 



a 2A- 2 / 3 C(4)[7V/C(3)] 4 / 3 and x c = 



a^A 1/3 C(4)[7V/C(3)] 4/3 for the axial and transverse di- 
rections (Fig. I2J, assuming an axial trap frequency 
u) z = Xuj±. a z = (Muj z ) and a± — (Mu>±) 
are the respective characteristic lengths for the harmonic 
trap along the axial and transversal directions. If we ap- 
proximate the effective area from the side view of a Bose 
gas as S w \J (z 2 )(x 2 ), we find the temperature depen- 
dence of S satisfies 



S = 



aS c (£) , T > T c , 



(7) 



with Sc = z'c x 'c ■ Once again, this is clearly consis- 
tent with earlier observations that the effective area is 
proportional to temperature above BEC and drops sud- 
denly below condensation temperature 0, l37l ] . 



III. SIZE OF AN INTERACTING BOSE GAS 

With atomic interactions, the Hamiltonian of our 
model system becomes 



H 



dr 



2 



in the second quantized form, where ^(^[^'''(r)] is the 
quantum field for annihilating (creating) an atom at lo- 
cation r and g — AirTi 2 a s /M denotes the interaction 
strength with a s being the s-wave scattering length. It 
is easy to show that the total number of atoms N is con- 
served since it commutes with the Hamiltonian. We then 
introduce a Lagrange multiplier fi to fix the total number 
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TABLE I: Parameters for various types of typical harmonic 
traps, in units of (27r)Hz. The cigar-shaped trap is for the 
experiment of Ref. Q. 





u x 


UJy 




Spherically symmetric 


65 


65 


65 


Cigar-shaped 


140 


140 


14 


Pancake-shaped 


30 


30 


300 



of atoms in our numerical calculations. In fact, [i is sim- 
ply the chemical potential. The target function we need 
to minimize then becomes the free energy 



G 



H 



dr 



UN 



2M 



V^ + Vext(r)-M * 



2 



(8) 



Adopting the standard mean field theory |23l I38j , we in- 
troduce a condensed (<£> = ($)) and a thermal component 
of a Bose gas through W = $ + <5\f'. After straightforward 
manipulations, we obtain the Gross-Pitaevskii equation 
for the two components as, 



d 

ih—$(f,t) 



ih—5$(r,t) = 
at 



TTTT V 2 + V cxt - n + g{n + n T ) 



2M 



~—V 2 + V cxt -fi + 2gn 



4>. 

(9) 

6V, (10) 



where n = nc + nr is the total density with nc = \&\ 2 
being the condensate density and n-j- the density of ther- 
mal atoms [given by Eq. (| 1 1 1) ] . In obtaining the above 
closed set of equations, we have used the Hartree-Fock 
(HF) approximation, which is both efficient for numeri- 
cal solutions and reasonably accurate for predicting the 
statistical properties of a Bose gas as confirmed by the 
excellent a gree ment between experiments and theoretical 

results [pii m m m m m . 

We will treat the thermal component through a semi- 
classical approximation, —iKS7 — > p, which gives its 
distribution in terms of a Bose-Einstein distribution in 
phase space {p, f } 



n T {r) 



dp 



1 



(2tt^) 3 e £ (p")/ fc BT _ i • 



(11) 



with e(p,r) = p 2 /(2M) + V cxt — fi + 2gn. Numerically 
we follow the standard procedure solving for the self- 
consistent solution of the coupled Eqs. JU}, CHJi an d 

CD MIH- 

We compute below the temperature dependent mean 
field ground state of an interacting 87 Rb atomic Bose gas 
inside a harmonic trap 

V ext (x,y,z) = -M{u>lx 2 +ujly 2 +ulz 2 ), (12) 



taking the atomic s-wave scattering length a s — 100. 4as 
with as the Bohr radius 39]. For spherically symmetric 
traps, we search for the ground state with different total 
number of atoms, N = 5 x 10 5 , 1 x 10 5 , and 1 x 10 4 . The 
temperature dependence of the effective width, defined 
as (r 2 ) = J drr 2 n(f), is shown in Fig. ^ Clearly a sud- 
den drop of the effective width for an atomic cloud occurs 
when temperature is lower than the critical temperature, 
as in the experiments Q, |3tJ • We further observe that 
the sudden decreasing of the reduced width near the crit- 
ical temperature has little dependence on the number of 
atoms. 

Figure ^ also reveals that the difference between an 
ideal gas and an interacting one increases with the num- 
ber of atoms, and with increasing repulsive interaction 
strength for T < Tq- An interesting feature we note 
is that the repulsive interaction causes the width of a 
Bose gas to increase at temperatures lower than the crit- 
ical temperature (T < Tc), yet it has little effect on the 
width at temperatures higher than the critical tempera- 
ture (T > Tq). The low temperature phenomenon is easy 
to understand in terms of a repulsive-interaction induced 
expansion of a Bose gas j23|. First, a condensate with re- 
pulsive interaction is larger in its size due to atom-atom 
repulsion; Second, the presence of a condensate pushes 
the thermal non-condensed cloud out, further increas- 
ing the width of a gas [3^, E3- At high temperatures 
(T > Tc) the effect of repulsive interaction becomes neg- 
ligible as the density of a Bose gas decreases dramatically 
with increasing temperatures. 

For cylindrically symmetric harmonic traps, we com- 
puted the ground states both for a cigar-shaped harmonic 
trap with an aspect ratio io z jio x = 0.1 and a pancake- 
shaped harmonic trap with lo z /lo x — 10. As before, we 
used TV = 5 x 10 5 . Of particular interest, we find excellent 
agreement with the experimental observations of Wang et 
al. Q for the above cigar shaped trap, at approximately 
the same transition temperature Tc ~ 230nK for their 
experimental parameters. In the experiment of Wang et 
al. 0| , the magnetic trap cannot be switched off without 
causing violent perturbations to trapped atoms. There- 
fore, in situ near resonance imaging [4l| was used instead 
of the conventional approach of shutting off the trap and 
letting a condensate expand [l5||. After condensation, a 
characteristic halo-like structures was observed in near 
resonantly diffracted light 0. Normally this cannot be 
taken as a convincing evidence for BEC. But for a cigar 
shaped trap as was used in their experiment, the halo 
is mainly along the transverse direction, while the axial 
length can still be determined with sufficient accuracy 
because it is relatively long, in fact much longer than 
the resonant optical wavelength; When the temperature 
dependence of the effective length for their experiment 
was plotted, a sudden decrease was observed, which was 
interpreted as an indication of the BEC phase transition 
jjl| . Our work thus provides a solid theoretical foundation 
for their experiment. As illustrated in the above figures, 
a sudden decrease of the effective length of the cloud was 
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observed when the temperature is lower than the tran- 
sition temperature, thus a corresponding decrease of the 
effective area. Similar to spherically symmetric traps, the 
difference between the effective widths of an interacting 
Bose gas and an ideal gas decreases when the tempera- 
ture increases and becomes negligible for T > Tc- We 
also note that the ratio between the transverse and the 
axial width remains about the same as the aspect ratio 
lj z /lj x for the interacting Bose gases. This is not sur- 
prising because this ratio approaches that of an ideal gas 
x cl z c — u z /u> x when T > Tc or that of the TF limit 
Xtf/ztf = w 2 /wj when T — > 0. 



IV. PARTITION OF ENERGY FOR AN 
INTERACTING BOSE GAS 

The total energy of a Bose gas can expressed in terms 
of its various partitions |23l |25| 



E = E { C K) + E c 



with 



E { r K) = [ dr— |V$(r)| 2 , 
c J 2M l Wl ' 

Ejp = J drV cxt (r)mr)\ 2 , 

E ( c ] = §/ dr\<S>(r)\\ 



E, 



(K) 



drdp p 2 /(2M) 



E, 



(V) 



(2ttK) 3 e<^f)/k B T _ 1 > 
drV cxt (r)n T (r), 

4 7) = 9 J drn 2 T (r), 
(I) - 2g f dr \<S>{r)\ 2 n T {r). 



E, 



CT 



E. 



(K) 



E^/ir,, and E^',^ denote respectively the kinetic 



(V) nrir i p (i) 

C/T' ^C/T> cUlu ^C/T 

energy, the trap potential energy, and the interaction 
energy partitions of the condensed/thermal components. 
E^Op denotes the interaction energy between the conden- 
sate and the thermal component. Figure [3] compares the 
temperature dependence for various energy partitions of 
an ideal Bose gas (blue lines) with an interacting 87 Rb 
gas (red lines) in a spherical harmonic trap. For an ideal 
gas, only kinetic and trap potential energy partitions are 
nonzero. In fact they are identically the same for ei- 
ther the condensate or the thermal component. The 
total energy of the thermal component exceeds that of 
the condensate when T > 0.3Tc- For an interacting 
gas, we see from Fig. that the kinetic energy of the 

condensate E^ is negligibly small though it still de- 
creases from a nonzero value at T — to zero at T — Tc 
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FIG. 3: (Color online) The temperature-dependent energy 
partitions of a Bose gas in a spherically symmetric harmonic 
trap at N = 5 x 10 s and uj — (2-7r)65 Hz. The upper panel 
illustrates the overall feature while the lower panel highlights 
the details of several small energy components. In the lower 
panel, the two blue dotted lines refer to the kinetic and (trap) 
potential energy for the thermal component of an ideal Bose 
gas. They are essentially identical. Similarly, the two solid 
blue lines denote the same energy components as for the con- 
densed part. The markers used for the various energy com- 
ponents of an interacting Bose gas (in red color) are denned 
in the legend. 



if we inspect the curve carefully. At low temperatures, 



E, 



3/2, which is a prediction from the Thomas- 
Fermi approximation. At high temperatures, E^ and 
E^P both approach zero as T — > Tc- Both kinetic and 

potential energies E^ and EjP of the thermal cloud in- 
crease rapidly as temperature increases (T < Tc), over- 
whelming the other energy partitions near approximately 
Tc/2, and increase linearly for T > T c . E { ^ ) is slightly 

larger than E^ because of the repulsive atomic interac- 
tion. They both approach their corresponding value for 
an ideal gas as T > Tc, as shown in Fig. [3J The inter- 
action energy of the thermal cloud E/p reaches its max- 
imum at ^< Tc and decreases slightly as T increases. 
We also observe from Fig. [31 that the interaction energy 
between the condensate and the thermal cloud Eqj* is 
largest near ~ Q.&Tc and decreases to zero at ~ Tc- 

Figure 01 illustrates the temperature dependence of the 
total release energy of a Bose gas from a spherical har- 



monic trap, £"r. = &c 



E. 



E: 



E: 



(I) 



see that essentially no difference exists between that of 
an ideal or an interacting gas. Comparing Figs. [3]with^ 
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FIG. 4: (Color online) The temperature-dependent release 
energy of a Bose gas for the same parameters as in Fig. 
The solid line is for an ideal Bose gas while circles denote an 
interacting 87 Rb one. 



atomic Bose gas. We found that the effective size can 
serve as a good indication for the presence of BEC. 
Our mean field theory gives several analytic results 
in the various temperature limits, and allows for a 
careful comparison between the temperature dependence 
of the various observable quantities between an ideal 
and an interacting Bose gas. In addition, the simple 
temperature dependence of the width of a trapped Bose 
gas can be potentially applied to infer the temperature 
of an atomic gas cloud |J, A detailed analysis 

of the energy partitions show that kinetic energy and 
trap potential energy of the thermal cloud represent 
the main contribution to the energy of the system at 
high temperatures, consequently they dominate the 
statistical properties of the system. Furthermore, our 
result provides an excellent theoretical support for the 
observation of atomic BEC in Ref. pj. 



we see the main feature of the release energy at high tem- 
peratures again follow the characteristic temperature de- 
pendence, i.e. E K oc (T/T c ) 4 if T < T c and E R oc T/T c 
if T > Tc, due mainly to the temperature dependence of 
the kinetic energy of the thermal component E^F ' . 

V. CONCLUSION 

In conclusion, we have investigated theoretically the 
temperature dependence of the average size of a trapped 
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